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P REFACE 


This  Memorandum  is  one  result  of  an  investigation  of  the  radar 
returns  from  missile  wakes.  Related  Memoranda  are  RM-3573-ARPA,  The 
Scattering  of  Electromagnetic  Waves  from  Plasma  Cylinders:  Part  I ,  a nd 
RM-3574-ARPA,  The  Scattering  of  Electromagnetic  Waves  from  Plasma  Cyl¬ 
inders;  Part  II. 

The  work  concludes  RAND's  study  of  midcourse  phenomenology,  con¬ 
ducted  under  the  Advanced  Research  Projects  Agency's  Defender  program. 
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SUMMARY 


This  Memorandum  investigates  the  steady-state  perturbations  or 
the  ion  and  electron  densities  caused  by  a  large-radius  conducting 
disk  moving  through  the  ionosphere  a*:  a  speed  intermediate  to  the  ion 
and  electron  thermal  speeds.  An  approximation  to  the  electron  density 
in  the  region  close  behind  the  body  is  found  by  solving  Poisson's 
equation,  neglecting  the  ion  density.  For  larger  distances  behind  the 
body  where  the  ion  and  electron  densities  are  almost  equal,  an  ex¬ 
pression  for  the  Fourier  transform  of  the  perturbed  ion  density  is 
obtained,  taking  account  of  the  ion  thermal  motion  and  scattering  in 
the  electric  anu  magnetic  fields,  but  neglecting  collisions.  The 
transform  is  inverted  to  obtain  an  approximate  expression  for  the 
space  dependence  of  the  perturbation  at  large  distances.  Finally, 
the  solutions  for  uifferent  body  shapes  and  boundary  conditions  are 
discussed. 
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I,  INTRODUCTION 

Let  us  consider  the  problem  of  the  steady-state  perturbations  in 
the  medium  caused  by  a  conducting  disk  moving  parallel  to  its  axis 

through  the  ionsophere  under  the  following  conditions: 

6 

1.  The  speed  of  the  bodv,  Vg,  is  about  10  cm/sec,  intermediate 

to  the  ion  thermal  snccd  (.V  10"*  cm/secl  and  the  electron  thermal 

7  1 

speed  (V>  ~  10  cm/sec). 

2.  The  radius  of  the  body,  R  ,  is  large  compared  to  the  electron 
Debye  length.  ~  1  cm). 

3.  All  ions  which  strike  the  body  disappear  as  ions. 

4.  The  Larmor  radius  of  the  ions  in  the  earth's  magnetic  field, 

Ps  ,  is  much  larger  than  the  body  radius  (based  on  the  speed  V  ,  R^  is 
of  the  order  of  50  m) . 

5.  The  ion  collision  frequency  in  the  medium  is  small  compared 
to  the  ion  cyclotron  frequency;  this  condition  holds  for  altitudes 
above  about.  130  km. 

6.  The  potential  of  the  body,  cp^,  is  not  much  larger  in  magni¬ 
tude  than  kT^  (T^  is  the  electron  temperature). 

Under  these  conditions ,  the  medium  can  be  xpected  to  be  quali¬ 
tatively  that  as  described  by  Gurevich. ^ ^  The  electrons  can  be  re¬ 
garded  as  being  in  equilibrium  in  the  potential;  because  of  the  great: 
speed  of  the  electrons  as  compared  to  that  of  the  body,  the  body 
appears  to  be  essentially  stationary. 

For  the  ions,  the  situation  is  much  different.  The  region  direct¬ 
ly  behind  the  body  will  be  essentially  evacuated  of  ions-  As  one  moves 
farther  back,  the  space  will  gradually  fill  up  with  ions  becau.se  o t 
their  random  thermal  motion  and  their  scattering  in  the  electric  field. 
Still  farther  back  (at  distances  larger  than  the  ion  Larmor  radius), 
the  ion  density  will  also  be  influenced  by  the  magnetic  field.  It 
should  be  pointed  out  at  this  time  that  while  tne  magnetic  field  is 
external,  the  electric  field  is  a  consequence  of  the  presence  of  the 
moving  body.  There  is  a  charge  on  the  body  itself  and  a  charge  den¬ 
sity  within  the  medium.  The  latter  results  from  the  different  influ¬ 
ence  of  the  moving  body  on  the  ions  and  electrons,  causing  different 
perturbations  in  the  electron  and  ion  densities. 
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The  problem  of  estimating  the  density  of  ions,  taking  account  of 
their  absorption  by  the  body,  their  thermal  motion,  and  their  scatter¬ 
ing  in  the  electric  and  magnetic  fie’ds,  wourd  appear  at  first  to  be 
tremendously  complicated.  However,  the  situation  is  somewhat  simpli¬ 
fied  by  the  fact  that  close  behind  the  body  where  the  perturbation  in 
the  ion  density  from  the  ambient  value  is  large,  electromagnetic  effects 
on  the  ion  density  are  minor,  and  the  ion  distribution  differs  little 
from  that  of  neutral  particles  having  the  same  mass,  temperature,  and 
surface  interactions  as  the  ions.  Farther  back,  where  the  effects  of 
the  electric  and  magnetic  fields  become  important,  the  perturbation 
in  the  ion  density  from  its  ambient  value  is  small.  Therefore,  in 
the  region  in  which  it  is  necessary  to  include  all  influences  on  the 
ion  density,  we  can  employ  small-perturbation  theory. 
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II.  GUREVICH  ZONES 


For  the  purpose  of  determining  the  electrostatic  potential, 
Gurevich^^  pointed  out  that  the  space  behind  the  body  can  be  regarded 
as  divided  into  zones. 

The  zone  close-.t  to  the  body  (denoted  by  Gurevich  as  the  maximum- 
rarefaction  zone)  is  one  in  which  the  ion  density  is  so  low  that  to  a 
first  approximation  it  can  be  neglected  in  the  determination  of  the 
potential. 

Farther  back  (at  distances  larger  than  about  R^( V^/V’  )  behind 

the  body),  the  ion  density  determines  the  potential;  except  for  terms 
2  2 

of  order  (assumed  here  to  be  small),  the  electron  density  equals 

the  ion  density. 

The  division  of  the  spac.  into  these  zones  can  be  seen  from 
Poisson:s  equation  for  the  potential  cp.  Using  dimensionless  coordi¬ 
nates,  with  the  body  radius  as  the  unit  of  length  and  the  origin  of 
the  coordinate  system  at  the  center  of  the  disk,  Poisson's  equation 
is  given  by 
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ne  '  ”l 
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where 
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1/P 


2 
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ecp/kT^,  a  measure  of  the  ratio  of  Che  potential  to 
the  kinetic  energy  of  the  electrons 
2  2 

X  /R  «  1,  i.e.,  the  ratio  of  Debye  length  to  body 

D  D 

radius  is  small 

the  electron  and  ion  densities,  respectively 

the  ambient  density  of  both  kinds  of  charged  particles 


The  conditions  on  $  are 


=  $  cn  the  body 
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0  at  infinity 
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3nd  $  aal  its  derivatives  are  continuous  (except  the  normal  derivative 
on  the  surface  of  the  body). 

Since  the  electrons  are  essentially  in  equilibrium  in  the  poten¬ 
tial,  the  electron  density  is  given  by 


n 

e 

n 

00 


$ 

e 


(2) 


The  ion  density  is  as  yet  to  be  determined.  However,  we  know 
that  there  will  be  a  region  close  behind  the  body  in  which  the  ion 
density  will  be  small  enough  so  that  rLj./n^  <  e’.  In  this  region,  as 
a  first  approximation,  we  can  neglect  n^/n^  altogether,  and  Poisson's 
equation  becomes 

=  e$  (3) 

PB 

Gurevich  pointed  out  that  in  the  maximum-rarefaction  zone,  which 

is  the  region  where  Eq.  (3)  is  applicable,  the  solution  ef  will  have 

2 

an  order  of  magnitude  of  1 / -  Denoting  $  within  the  maximum-rare¬ 
faction  zone  as  the  zero-order  approximation  to  §  as  given 

rlK4i  rlKZ 

by  Gurevich  is 


(4) 


Gurevich  did  not  obtain  the  next  approximation  in  this  region. 

Outside  the  maximum-rarefaction  zone,  where  n^/n^  becomes  com- 

A 

parable  to  e’ ,  Gurevich  pointed  out  that  a  solution  to  Poisson's  equa¬ 
tion  can  be  found  by  an  iteration  procedure.  This  is  effected  by 
writing  Poisson's  equation  in  the  form 


=  In 


nl  ,  1  2 

7T  +  ~V  * 

.  -  pB  - 


(3) 
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For  the  zero-order  approximation  to  $,  the  ion  and  electron  densi¬ 
ties  are  equated: 


-  in 


n 


The  next  approximation  is 


(o) 


(7) 


and  so  forth.  At  large  enough  distances  behind  the  body,  the  zero- 
order  approximation  (Eq.  (6))  suffices.  The  problem  in  this  region, 
then,  is  to  determine  n^..  Gurevich  did  this  by  first  treating  the 
ions  as  neutral  particles  and  then  taking  account  of  the  magnetic 
field.  He  did  not  consider  the  effects  of  the  electric  field  on  the 
ion  density. 

In  this  Memorandum  an  attempt  will  first  be  made  to  find  a  solu¬ 
tion  for  ky  considering  a  highly  idealized  model  in  which  the 

ions  have  zero  temperature  and  in  which  the  influence  of  the  electric 
and  magnetic  fields  on  the  ion  motion  is  ignored.  The  solution  for 
this  model  should  approximate  the  actual  solution  within  the  maximum- 
rarefac^ion  zone. 

Second,  an  attempt  will  be  made  to  estimate  the  ion  density,  tak¬ 
ing  account  of  the  thermal  motion  of  the  ions  but  not  their  scattering 
in  the  electric  and  magnetic  fields-  Since  this  approximation  to  the 
ion  density  is  a  fairly  good  one  for  distances  behind  the  body  less 

than  (V  /V )R  ,  it  can  be  used  to  define  the  "boundary”  separating 
BIB 

the  maximum -rarefaction  zone  from  the  rest  of  the  space  behind  the 
body.  This  boundary,  which  is,  of  course,  fictitious  and  merely  a 
mathematical  convenience,  will  be  defined  as  the  surface  upon  which 

$MR? 

the  ion  density  is  equal  to  e  • 

Finally,  the  ion  density  will  be  considered,  taking  account  of 
the  electric  and  magnetic  fields,  in  the  region  where  the  ion  and 
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electron  densities  are  essentially  equal  and  Eq.  (6)  is  applicable. 

The  effect  of  collisions  will  be  neglected;  it  is  understood,  there¬ 
fore,  that  our  solution  is  invalid  for  distances  behind  the  body  larger 
than  the  collision  mean  free  path. 
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III.  SOLUTION  WHERE  TEMPERATURE  AND  FIELD  SCATTERING 
OF  IONS  ARE  NEGLECTED 


If  the  thermal  speed  and  scattering  of  the  ions  in  the  electric 
and  magnetic  fields  are  neglected,  then  the  ion  density  will  be  zero 
ir.  a  semi-infinite  cylinder  of  radius  RD  behind  the  body  and  equal  to 

D 

the  ambient  density  everywhere  else.  Using  cylindrical  coordinates 
Z,  p,  and  0,  Eq.  (1)  becomes 


for  p  <  1 


(8) 


for  p  >  1 


Negative  Z  refers  to  the  region  behind  the  body. 

The  boundary  conditions  are  the  same  as  those  for  Eq.  (1)  except 
for  the  condition  on  $  as  Z  -*  -®.  An  assumption  that  $  -*  0  as  Z  —  -co 
implies  that  the  electron  density  becomes  equal  to  its  ambient  value 
as  Z  -*  -oo.  However,  for  identically  zero,  the  ion  density  is  zero 

in  the  entire  semi-infinite  cylindrical  region  Z  0,  p  <  1.  A  large 
charge  separation  such  as  this  is  physically  unreasonable.  We  would 
expect,  rather,  that  the  electron  density  would  remain  depressed  below 
its  ambient  value  throughout  the  ion  ''hole,"  implying  a  nonvanishing 
negative  potential  behind  the  body  extending  to  infinity  for  p  <  1. 
Indeed,  we  would  expect  that  far  behind  the  body,  $  becomes  essen¬ 
tially  independent  of  Z.  Therefore,  the  condition  that  $  -*  0  as 
Z  -*  -oo  will  be  replaced  by  the  physically  sensible  condition  that  as 
Z  -*  -“,  $  becomes  a  function  of  p  only. 

Therefore,  we  will  assume  that  $  can  be  written  in  the  form 


f  -  -u(p)  ~-Sl|q  +  w(p,  Z) 


(9) 
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where  w(p,  Z)  ->  0  as  jzj  -»  ®,  and  where  u(p)  satisfies  the  following 
equations  and  boundary  conditions: 


lid/  du\  -u 
2  p  dp  \p  dp  J*  "e 

PB 


P  <  1 


1_  i  ±_  (p 

2  p  dp  l  P  dp  i 

PB 
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-i; 
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p  >  1 


u  -*  0  as  p  -»  33 


(10a) 


(10b) 


and  u  and  du/dp  are  continuous  at  p  =  1.  The  function  w  will  satisfy 


r  ? 
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L  3”  i 

j)z2  P 
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\p  ^  /_ 

Z  >  0,  all  p  (11a) 


*u,  w  i  •. 
e  (e  -  1) 


Z  <  0,  all  p  (Hb) 


(w  -  u) 


Z=0- 


wz=o+ 


p  <  1 
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(u  -  u)z=o-  =  uz=o+ 


) .. 


OW  \ 

\  'z/z=o- 


dw 

dZ, 


z=of 


We  will  first  look  at  the  function  u(p),  which  determines  the  electron 
density  far  behind  the  body,  and  then  examine  the  function  w. 


SOLUTION  FAR  BEHIND  THE  BODY 


Interior  Solution  (Arbitrary  pD) 

*  D 


Equation  (10a)  can  be  solved  analytically: 


-u 

e 


8 _ a _ 

2  /  .  2\ 2 
PB  l  1  *  <*p  I 


P  -  1 


(12) 


where 


and  u  is  the  value  of  u  on  the  axis  (p  «  0).  The  value  of  u  on  the 
o 

axis  is  as  yet  unknown;  it  is  determined  by  matching  the  inside  (p  £  1) 

solution  to  the  outside  (p  >  1)  solution  at  p  -  1.  In  terms  of  a  and 

p  ,  u  and  its  derivative  at  p  =  1  are  given  by 
B 


i 


-U, 


=  2pB  ore 


(13a) 


(13b) 
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whcre  u^  Is  the  value  of  u  at  p  =  1. 

We  will  use  Eqs.  (13a)  and  (13b)  later  on.  One  thing  we  can  see 
immediately  from  the  i nside  solution  is  that  a  must  be  less  than  one 

in  order  for  e  U  to  be  finite  for  all  p  £  i.  Thus,  n  /n  must  be  less 

2  s  e  “ 
than  8/pp  on  the  axis.  As  woulch  Jbe  expected,  the  electron  density 

near  the  axis  gets  very  small  as  the  radius  of  the  "hole"  gets  very 

large. 


Potential  at  p  =  1 

Before  attempting  to  find  a  solution  for  p  >  1,  let  us  first  see 
how  much  information  can  be  obtained  from  Eqs.  (10b)  and  (13)  without 
getting  an  explicit  solution  co  Eq.  (10b).  Let  us  rewrite  Eq.  (10b) 
as  follows: 


l_  d/jj 

2  ,  2 
PB  ^ 


-u 


1  du 

P  dp 
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(14) 


Multiplying  both  sides  of  the  equation  by  du  and  integrating  from  0 
to  u^ ,  we  get 


r- 
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1  /du 
■B 
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jl  r 1  _i _ du 

2  J  p(u)  dp 
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(15) 


Combining  Eqs.  (13b)  and  (15),  we  get 


-u , 


u 


(1  -  r 


1 


p2  0  p(u)  dp 


du  , 
T-  du 


(16) 


Since  du/dp  is  negative,  the  right-hand  side  of  Eq.  (16)  is  positive. 
This  means  that  since  a  <  1,  u,  must  be  less  than  on.  for  all  values 

1  2  _ui 

of  p  .  Thus,  for  large  Pg,  the  term  p0  e  1  will  be  large,  and  Eq. 
(13a)  becomes 
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,  V2 

4  e _ 
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(17) 


Moreover,  since  both  ( 1  -  cv)  and  the  right-hand  side  of  Eq.  (16)  are  of 
rder  1 / p  ,  u  differs  from  unity  by  a  term  of  order  1 / p  for  large  p  . 

U  I  D  D 

Thus,  we  have 
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Interior  Solution  for  Large  p 

D 


Knowing  a,  Eq.  (12)  can  now  be  rewritten  for  large  po: 
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-u 
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(19) 


Since  the  electron  density  far  behind  the  body  is  given  by  n^e  u, 
it  can  be  seen  from  Eq.  (19)  that  / n^  is  very  small  in  the  interior 

of  the  large-radius  ion  "hole."  On  the  axis,  n  /n  is  equal  to  8/p^ 

e  as  o 

and  rises  fairly  slowly  with  increasing  radius  up  to  a  distance  of  a 
few  Debye  lengths  from  the  edge  of  the  cylinder.  Thereafter,  n  /n 
rises  rapidly  until  p  *  1.  At  p  =  1,  it  is  equal  to  about  1/e,  pro¬ 
vided  that  p_  is  large.  Thus,  a  simplified  picture  of  the  region 

D 

p  5  1  far  behind  the  large-radius  body  is  one  of  a  cylindrical  hole 
devoid  of  ions  and  low  in  electron  density,  except  for  a  sheath  of 
thickness  of  a  few  Debye  lengths  at  the  edge.  The  electron  density 
in  the  sheath  is  given  by 


n 

e 


n 

ac 


(20) 


where  p  =  (p  -  l)p  .  Thus,  the  electron  density  inside  the  cylinder 

D 

near  the  edge  depends  upon  the  distance  jpj  in  De^/e  lengths  from  the 
edge  of  the  cylinder,  not  on  the  value  of  for  large  . 


Solution  for  p  ^  1  for  Large  o 

D 


For  large  p  ,  it  is  convenient  to  rewrite  Eq.  (10b)  as  follows: 

D 


U"(p)  + 


PB^ 


- r-  U  '  (  0)  =  1 

P/P  J 


-U 


(21) 


Neglecting  terms  of  order  1 / p  ,  Eq.  (21)  becomes 

D 


u"(p)  =  1 


-u 

e 


V.  —  “/ 


Ore  integration  gives 


(23) 


so  that 


u, 


1_  j 

V?  J 


ou 


/  -u 

u  .  /  e  +  u 
V 


*  p 


(24) 


Since  u,  ~  1  and  u  is  a  monotonical ly  decreasing  function  of  p,  we  can 

^  vi  3 

expand  the  exponential  e  .  Keeping  terms  up  to  u  ,  we  get 


u  ~ 


(1 


(25) 


Since  u  decreases  exponentially  with  p,  the  electron  density  becomes 
essentially  equal  to  its  ambient  density  within  a  few  Debye  lengths 
from  the  edge  of  the  cylinder.  Moreover,  we  see  that  the  electron 
density  near  the  edge  outside  as  well  as  inside  the  cylindrical  hoi'’ 
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depends  upon  the  distance  in  Debye  lengths  from  the  edge  and  not  on 
the  value  of  p  ,  as  long  as  p  is  large. 

D  D 

The  electron  density  and  tue  absolute  value  of  the  charge  density 
are  plotted  as  functions  of  p  in  Figs.  1  and  2  for  p„  =  20,  50,  and 

O 

100.  Figure  2  shows  that  the  absolute  value  of  the  cnarge  density 
(divided  by  en  )  becomes  small  within  a  few  Debye  lengths  from  the 
edge  of  the  cylinder,  as  would  be  expected. 

ELECTRON  DENSITY  NEAR  THE  DISK 

In  order  to  obtain  the  potential  near  the  body,  we  need  to  esti¬ 
mate  the  function  w(p,  Z) . 

Solution  Inside  the  Hole  for  Large  pD 

iJ 

We  will  find  an  approxi.nate  solution  to  Eq.  (lib)  inside  the  hole 
for  large  by  first  rewriting  the  equation  in  terms  of  p  and  a  new 

D 

variable,  C,  where 


(u  -u)/2 

C  =  -Ze 


(26) 


Equation  (lib)  becomes 


where 


|  ■—  =  eW  -  gx(p,  0  -  g2(p,  C) 


(27) 


8i  “ 


1  bw 
9  C  xr  + 


1  2  b  9 


u-u  /  0 
l  +  1  (i-S  +  1 

•  W  > 9p 


16 


?  (u-u  )/2 

Sw  b  'w  ,  e 

SC  Sc2  p 


b  ( bw^ 


the  charge  density  far  behind 
zero  ion  temperature 
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The  function  g,.  vanishes  on  the  axis  for  ail  '  and  on  the  surface  of 
£ 

the  disk  for  all  p  >  1.  For  and  p  small  but.  finite,  can  be  ex¬ 

pected  to  be  small  in  magnitude  The  function  g^  will  be  of  order 
unity  for  large  '  when  w  becomes  small.  For  small  f,  $  =  $  (in¬ 
dependent  of  p) ,  so  that 


1  o  /  dw\  1  d 

-  —  1  p  ~  u  — 

P  Bp  \  dp  /  p  dp 


u 


Thus 


as  f 


Therefore s  for  small  Q,  both  and  g^  will  be  small  in  magnitude,  be- 

w 

coming  identically  zero  on  the  surface  of  the  disk.  The  term  e  ,  on 

the  other  hand,  will  be  large  for  small  f  and  p  if  p  is  large  and  if 

o 

is  of  order  unity.  If  we  restrict  ourselves  to  values  of  -$  £  1, 

b  d 

then  e  >  1  for  all  p  ^  1.  We  also  know  that  for  large  £,  w  will 
approach  zero. 

Thus,  it  seems  reasonable  to  replace  Eq.  (27)  by  the  following: 


i  a2 
1  d  w 


w 

e 


£  *  £*(p)  (28) 


W  =  Ig  +  u ( p)  =  Wg(p)  £  -  o 


w 


dw 

d£ 


0 


£  >  C*(p) 


The  function  r  (p)  is  as  yet  unknown;  it  is  the  surface  upon  which 
both  w  and  dw/d'  vanish  and  is  determined  from  the  solution  to  Eq.  (28). 
One  integration  of  Eq.  (28)  yields 
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7 

1  dw 
16  d 


(29; 


The  constant  of  integration  Was  chosen  to  cause  dw/'d '  to  vanish  when 
w  =  0.  One  more  integration  gives 


1  = 


v<vT 


tan  2C, 


L+yfTT 


tan  2' 


(30) 


or 


w 

e  = 


WB  2 
e  sec  2u 


(1+^v7 


tan  2£ 


r  <  1  -1  _WB/2  * 

Q  £  -  cos  e  -  b 


The  function  £  (p)  is  then  given  by 


*  1  -l(  (u‘V/2) 

C  (P)  =  o  cos  \e  / 


(31) 


Written  in  terns  of  Z  rather  than  £,  Eq.  (31)  becomes 


*  *  (u-u  )/2 

Z  (p)  -  -C  e 


z*(p)  .  .  Li. °£2J.  cos-f^ 


•V2 


PB  (1  -  Q'p2) , 


(32) 


with  a  given  by  Eq.  (18). 

'fc 

The  value  of  -Z  will  vary  from  approximately  tt/4  at  p  =  C  to  a 

-  1.  We  can  see  _rcm  Eq.  (32)  that  the  arc  cosine 


value  near  2ero  as  p 
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2 

term  will  be  about  tt/2  except  for  p  near  1;  for  o  near  1,  the  1  -  yp 

factor  will  become  very  small,  since  O'  is  near  unity.  Therefore,  an 

excellent  approximatio.  of  Z  (p)  to  Eq.  (32)  for  large  p  and  -$ 

approx  15  d 

of  order  unity  is 

-Z*  (p)  a:  7  (1  -  p2)  (33) 

approx  4  v  ' 

independent  of  p  and  $  .  As  an  illustration,  both  Eq.  (32)  and  Eq. 

D  D 

(33)  are  plotted  in  Fig.  3  for  p  =  100  and  -$  =  1. 

D  15 

w 

We  will  not  attempt  to  obtain  a  more  accurate  solution  for  e 
than  that  given  by  Eq.  (30).  The  error  involved  in  replacing  Eq.  (27) 
by  Eq.  (28)  becomes  large  when  w  becomes  small;  however,  this  error  is 
not  likely  to  be  larger  than  other  errors  inherent  in  our  approximate 
theory.  Furthermore,  it  should  be  noted  that  the  error  in  our  approxi¬ 
mations  is  only  significant  if  we  are  interested  in  derivatives  of  the 
electron  density.  The  electron  density  itself  is  proportional  to  eW, 
and  when  w  is  small,  errors  in  w  of  large  percentage  will  not  matter 
very  much. 

Thus,  we  have  for  the  potential  behind  the  disk  inside  the  large- 
radius  cylindrical  hole 


e  =  e 


G  s  G 


(34) 


*MRZ  -u 
e  3  e 


G  *  G 


In  terms  of  Z  and  p,  this  gives,  for  p  not  too  near  1, 


-20- 


rMRZ 


'B 


Ulp,  V2 

1  +7r-e 


i2  ys 


ZIPB  V2<1 

e 


$MRZ  -u  2  n  _  -u  l(  2Z  \ 
e  e  esc  2r  =  e  esc  i  - -  | 

V  1  -  P2  / 


-$o/ 2 


VT  -(pr 

e  <  |Z|  <  | Z  ( p) 
h  B 


MRZ  -u 
e  =  e 


|z|  *  |z*(P)| 


(35) 


Figure  4  plots  Eq.  (34)  along  the  axis  p  =  0  for  p_  =  100  and  =  -1. 

b  D 


Solution  Outside  the  Hole  for  Large  p 
_ , _ _  _ _  B 

For  Z  >  0,  w  *  $.  If  is  not  larger  than  one,  then  jw|  will 
be  less  than  one  for  the  entire  region  Z  >  0. 

For  Z  <  0,  w  =  $  +  u.  We  have  already  determined  that  u  <  1  out¬ 
side  the  hole.  If  -$  ~  1,  |w|  will  also  be  small. 

Therefore,  we  can  linearize  Eq.  (1)  outside  the  hole  for  both 
positive  and  negative  Z: 


l\2 

0  w 

1 

3 

(  dw\ 

-72 

oZ 

m. 

+  — 

P 

dp 

\P 

-  w  = 

Z  >  0 


1 


P 


2 

B 


d2w  1  d 
3z2  P 


-u 

e 


w  =  0 


Z  <  0 
P  >  1 


(36a) 


(36b) 


The  function  e  U  in  Eq.  (36b)  is  equal  to  1/e  at  p  =  1  and  becomes 
equal  to  one  within  about  a  Debye  length  from  the  edge  of  the  cylinder. 
The  solutions  to  Eq.  (36)  that  satisfy  the  boundary  conditions  that 
jwj  vanishes  as  p  and  |z|  become  infinite  decrease  exponentially  with 


3-/9  7T/6 

iZi  (rodians) 


/4 . 5  7TV 


ectron  density  behind  the  body  at 
p  =100  Debye  lengths  ( <E>  *  -1) 

B  B 
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the  coordinates  measured  in  Debye  lengths  from  the  boundaries  of  the 
semi-infinite  cylindrical  hole.  In  light  of  the  approximations  in¬ 
volved,  finding  exact  solutions  seems  pointless.  It  suffices  to  say 
that  within  one  or  two  Debye  lengths  from  the  bouidaries  of  the  cylin¬ 
drical  hole,  the  electron  density  is  equal  to  its  ambient  value. 
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IV.  BOUNDARY  Z_(p)  BETWEEN  THE  GUREVICH  ZONES 

G  _  _ _ _ _ 


Between  the  maximum- rarefaction  zone,  in  which  the  zero- ion- tem¬ 
perature  solution  for  the  electron  density  should  be  a  good  approxima¬ 
tion  to  the  actual  solution,  and  the  region  in  which  equating  the 
electron  and  ion  densities  provides  a  good  approximation  to  the  actual 

solution,  there  exists  a  transition  region  in  which  neither  approxima- 

* 

fcion  is  valid.  Within  this  transi  tion  region,  the  terms  e\,  /n  , 

2  2  x 
and  (l/p^)v  $  are  all  of  the  same  order  of  magnitude,  and  we  cannot 

justifiably  neglect  any  one  of  the  three  terms. 

The  surface  Z  =  2  (p),  upon  which  the  zero-ion-ten nerature  solu- 
$  G 

tion  fur  e  is  equal  to  the  ion  density  (divided  by  n  ) ,  will  lie  in 

Oft 

the  heart  of  this  transition  region.  Due  to  the  imprecise  nature  of 
this  "boundary,"  we  will  find  Z^(p)  using  the  simplified  approximation 
to  (which  we  shall  call  n ynj  obtained  by  treating  the  ions  as 

neutral  particles.  The  electric  field  will  not  alter  the  order  of 
magnitude  of  the  ion  density;  including  its  effect  on  for  the  pur¬ 
pose  of  obtaining  Z_(p)  would  simply  shift  the  surface  Z  =  Z  (p) 
slightly  within  the  transition  region.  The  magnetic  field  at  the  dis¬ 
tances  we  are  considering  here  (lz  ]  <<:  P  /R  )  is  irrelevant. 

If  the  scattering  in  the  electric  and  magnetic  fields  is  neglect¬ 
ed,  the  ion  density,  based  on  a  Maxwellian  distribution  at  infinity, 
is  given  by 


(37) 

where 


Z  -  -Z  TT  Z  <  0 

B 

and  where  I  is  the  hyperbolic  Bessel  function  finite  at  the  origin. 
For  p  =  0,  Eq.  (37)  becomes 


(p  -  0)  =  /  exp  - 


~2  :.2  rp — pr~i  yh  j  e  ' 

Z  +VI  /VBJWZ  -VVB  V  *1 


where 


■yi  =  V 


V  Z2  +  v2/v2 


—  (p  -  0) 


z  +  vr/v; 

1  15 


z2  +  v- /V2 


•z  >  vj/v2 


The  boundary  Z_,(p)  between  the  Gurevich  zones,  as  we  define  it, 
is  found  by  equating  nj /rv  as  given  by  Eq.  (37)  to  e  as  given  by 
Eq-  (19).  In  Fig.  5,  ^(p)  is  plotted  for  =  100  and  V2/V2  =  0.04, 
0-01,  and  0.  The  latter  case  corresponds  to  letting  VT/VZ  ->  0  m 

^  Id 

Eq.  (37)  while  keeping  Z  finite.  If  we  do  this,  Eq.  (37)  reduces  to 
the  much  simpler  expression 


^  -  2e-p2/Z~2 


f°  xe-*2  I  (H  x)  dx 

1/2  °\  z  / 


We  can  see  from  Fig.  5  that  for  the  parameters  we  have  chosen,  the 
three  curves  are  close  together.  In  general,  we  can  establish  a  cri¬ 
terion  for  replacing  Eq.  (37)  by  Eq.  (39)  for  the  purpose  of  finding 

Z  (p),  by  looking  at  the  expression  for  Z  (0)  : 

G  G 


V°>  / _ 1 _ \  _  "uo  _  _8 

t  9 — |  ^9  9  o  j ®  “  9 

ZG<°)+VVB  \ZG(0)  +VX/Vb/  PB 


2  2  ?. 

If  V  /V  is  large  compared  to  ;n  (p'/8),  then 

D  I  D 
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Z^(0)  «s( 


,1  - 


2  V 


fB 

2  /n  8 
B 


2\ 


) 


Jin 


W* ’> 


(41) 


will  not  depend  much  upon  V^./V  . 

We  recall  that  the  maximum- rarefaction  zone  is  divided  into  two 
subzones:  |zj  £  ]z  (p)|  and  jzj  t>  |z  (p)L  Therefore,  our  lowest 
approximation  to  $  for  Z  <  0  is  given  by 


=  -u(p)  +  w(p,  z) 


#n  "  -j<p> 


hi  *  iz*(p)i 


* 

Z  (P) 


P  <  1 


(42) 


Vg  ~ 

Z !  ^  —  ZG(p)  P<1 


"i 

§  =  £n  — 

III  n 

OD 


all  other  negative  Z 


The  location  of  the  three  regions,  I,  II,  and  III,  is  shown  in  Fig.  6 

for  pg  =  100  and  V^/V^  -  10-  Figure  7  shc-?s  few  equipotentials  baseci 

on  our  zero-order  approximations  to  $  for  p  =  100,  V  /V_  =10,  $  *  -1, 

.  ,  B  B  I  B 

and  |Z|  <  VB/V  . 


Fig. 7  —  Equipotentials  for  p&  =  100,  —=  10,  and<£0s-l 
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V.  ION  DENSITY  FOR  2  >  1 


In  order  to  determine  the  ion  density,  let  us  look  at  the  collision- 
free  steady-state  Boltzmann  equation  for  the  ion  distribution  function 

fj  (S,  V): 


(V  -  VB)  • 


2 

T  V, 
e  I 

T  2 
I 


7 9 


Vi 


(V  x  e  ) 
o 


S 


(A3) 


where  e  =  B/B  is  the  unit  vector  in  the  direction  of  the  magnetic 
o 

field  B,  and  p  =  (V  /eB/nijC)  1/Rfi  is  the  dimensionless  ion  Larmor  ra¬ 
dius  based  on  the  speed  VD.  The  right-hand  side  of  Eq.  (43)  (repre- 

D 

sented  by  S)  is  a  source  term  representing  the  change  in  the  distribu 
tion  function  due  t£  the  uptake  of  ions  by  the  disk;  S  can  be  written 
as 


[1  +  sign  (1  -  p)l  [1  +  sign  (V  -  V  )] 

5  =  5(Z)(VZ  -  VB)  - - - 2 - ~ -  fI  (44) 

(2) 

We  will  use  the  usual  technique  of  Fourier  transforms  to  obtain 
a  solution  of  Eq.  (43)  in  the  region  where  small-perturbation  theory  is 
applicable  ( j n^  -  n  | /n  <  1).  The  result  for  the  Fourier  transform  nR 
of  (nT  -  n  )/n  ,  based  on  a  Maxwellian  distribution  at  infinity,  is  given 

N  T  co'  <x 


\  = 


T 


e 


$K(1  +  iN) 


2n 


J^H) 


(45) 


where 
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<f  =  the  Fourier  transform  of  i 
K 

K  =  the  Z -component  of  K 


K 


and  where 


N  = 


d1]  exp  - 


2  n2 

^11  ^  2 

2  +  X1 


1  -  cos 


JL 


?lkz 


iTl 


= 


kii  vi 


K  = 


V2  2  w2 

kipl  !i 
2  2 


— ft 

(K 


2 


We  can  relate  ^  t  n^  by  means  of  the  Fourier  transform  or  Poisson’s 
equation.  However,  we  are  here  concerned  with  distances  that  will 
allow  us  to  equate  the  electron  and  ion  densities.  Therefore,  we  have 


*K 


°K 


-2tt 


JX(H) 

~K~ 


N 


1  +  -f-  (1  +  iN) 
I 


(46) 


Before  considering  the  inversion  of  the  Fourier  transform,  it  would 
be  useful  to  consider  ry  itself  under  various  circumstances,  since 
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for  some  problems  we  are  more  interested  in  than  in  n  -  n^.  (For 
example,  in  the  region  where  small-perturbation  theory  is  applicable, 
the  scattering  cross  section  of  radiation  from  the  wake  is  proportional 
to  ,  where  the  vector  K  refers  to  the  difference  in  the  propaga¬ 

tion  vectors  of  the  scattered  wave  and  the  incident  wave.) 

FOURIER  TRANSFORM  OF  PERTURBED  PARTICLE  DENSITY 

We  will  investigate  the  integral  N  in  various  limits. 


PL  -  00 

In  the  limit  of  infinite  Larmor  radius  (negligible  magnetic,  field), 
we  can  expand  the  cosine  in  the  exponent  of  the  integrand: 
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fTZjl 


expjx2  -  )  dx 
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(47) 


The  second  term  can  be  approximated  for  p,  small  and  large: 
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Thtis ,  when  the  magnetic  field  is  negligible,  is  given  by 


-2n 


J1(H)  j 
h  K. 


V 

Kz  ^  v'2 

VB 


°K 


1<H)  /T, 


i  v _ B  1  „2  ^  VI  2 

“  2  k  KV  1  +  T  /I  KZ  2  * 

1  01  VB 


If  the  electric  field  had  been  neglected,  Eq.  (49a)  would  be 

affected.  The  effect  of  the  electric  field  for  <  V2/\2  m2  (Eq. 

Z  X  B 

is  to  reduce  |n^|  by  a  factor  1/(1  +  T  /T_) . 

2  «  1;  K  •  e  4  0 

_ o _ 

2 

For  small  values  of  X  ,  the  integral  N  is  approximately 


N 


s=a 
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so  that 
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For  the  special  case  of  identically  zero  y,K  =  K.||),  we  have 


M- 


2 

I! 


2<%  •  V 


(49a) 


(49b) 
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so  that 
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(53b) 


K  .  e  -  0;  e  x  VD  *  0 
_ o _  o  B 

If  K  •  e  -*  0.  and  V„  is  not  parallel  to  B,  then  m- i i  =  0-  For  this 
o  B  II 

case,  we  write 
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so  that 
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(54) 


(55) 


Whenever  p  K  is  an  integer,  N  becomes  infinite.  However,  the  integral 

L  Li 

N  appears  in  both  the  numerator  and  denominator  of  Eq.  (46),  so  that 
n^  itself  is  finite  and  is  given  by 


CI  J1(H) 


*K  ■  2^L  T  Hi  PL 


KZPL  =  * 


(56) 


where  i  is  an  integer. 
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For  plK 
Eq.  (51a)* 

V_  Parallel  to  B;  k"  <  1 

D  _ 

If  the  body  is  moving  parallel  to  the  field,  p,||  is  independent 
of  K  and  is  equal  to 


2 

z  not  equal  to  an  integer ,  and  \i  small,  n^  is  given  by 


M.?.  -  (V^/2vJ)  «  1 


(57) 


Equation  (51a)  will  be  applicable,  with  n  <  l/pT,  V  /V  «  1,  so  that 

i-»  J.  B 


\  **  K, 


■2ni  J1(H) 


H 


-ni 


(58) 


For  motion  parallel  to  the  magnetic  field,  n^  will  go  to  infinity  as 
1/KZ  as  Kz  -  0. 

SPACE  DEPENDENCE  OF  PERTURBED  ION  (OR  ELECTRON)  DENSITY  FOR  Z  >  1 

The  major  difficulty  in  inverting  the  Fourier  transform  is  pro¬ 
vided  by  the  factor  1  / [  1  +  (1  +  iN)Te/Tj.  j,  which  cs  a  consequence  of 
the  electric  field..  However,  we  have  seen  that  for  large  Z  (small 
^Z^B^l))  ’  t^e  ^nfe8ral  N  is  small  in  absolute  value,  and  the  effect 
of  the  electric  field  is  to  reduce  the  magnitude  of  the  perturbation 
by  a  factor  of  approximately  1/(1  +  Tf/TJ.)  (except  for  motion  parallel 
to  the  field,  which  we  will  treat  separately). 

Therefore,  we  will  approximate  n^  for  large  values  of  2T  by 


•2n 


V*) 


N 


1  +  T  It 
e  I 


(59) 


The  inversion  for  p  «  p  and  p  «  Z  can  be  performed  in  a  straight¬ 
forward  manner.  The  result  it  given  by 
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ril 


n 

CO 


n 


as 


(1  +  Te/TI>(2pLvJ/vJ)|sin  Z/2ql\Jz2  sin2  *  +  4p£ 


2  2 

cos'  &  sin  (Z/2p, ) 

ij 


(60) 


where 


M 

2’l 


and 


COS  'll 


^  1 


For  jzj/2p.  <  1,  Eq.  (60)  becomes 

li 


00 


-~_-I _ l _ 1 _ 

1  +  T  /T  1  +  T  /T 

e  I  Z  el 


n 


n 

TO 


(61) 


As  expected,  for  distances  smaller  than  the  ion  Larmor  radius,  the  mag¬ 
netic  field  does  not  affect  the  solution.  Thus,  for  (V_/VT)R  <  | 3 |  <  R  , 
where  j z (  is  the  actual  distance  behind  the  body,  we  have  an  inverse- 
square  potential. 

For  |z|/2p^  >  cot  ’lr,  we  have 


"i 


-  n 


n 


as 


(1  +  T  /T  ) 
el 


sin  Z/2pJ(2pL  V^/VgjZ 


(62) 


sin  >if 


For  distances  large  compared  to  R,  ,  the  potential  oscillates  about  a 
mean  which  decays  as  1/jzj. 
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For  *  *  0,  F.q.  (p9)  does  not  apply  for  |zj  >  p  .  Tor  motion 

*>  ^ 

parallel  to  the  field,  the  parameter  is  independent  of  h  and  is 
small  for  all  K^.  Therefore,  we  ha  re  seen  (1  Jr  iN)  -*  0,  and  we  have 


J,(H) 

1L  ~  -2^  — —  -  -  $ 

1C  ~  hK„ 


(63) 


-  n 


1 


n 


(4pl  vi/vb}  sin2  <z/2pi/ 


lz!  >  pl 


(64) 


The  oscillations  in  $  for  ir  =  0  and  for  jzj  >  p  are  about  a  nonzero 
mean  which  is  independent  of  Z,  implying  an  extension  of  the  perturba¬ 
tion  to  infinity.  In  practice,  the  magnetic  field  is  never  strictly 
uniform,  and  maintain!  ng  a  value  of  identically  zero  is  impossible. 
Moreover,  the  effect  of  collisions  (which  we  have  completely  neglected) 
would  be  to  restore  the  medium  to  its  unperturbed  state  at  some  dis¬ 
tance  behind  the  bcdy. 
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VI.  EFFECT  OF  BODY  SHAPE  AND  MATERIAL 


Throughout  this  Memorandum  it  has  been  assumed  that  the  body  is  a 
conducting  disk.  It  is  clear,  however,  that  for  distances  benind  the 
body  larger  than  the  body  radius,  we  would  obtain  essentially  the 
same  solutions  for  any  body  with  a  circular  cross  section  in  the  plane 

—4 

normal  to  V  . 

D 

Moreover,  the  solutions  behind  the  body  for  jzj  ^  1  do  not  depend 

upon  the  potential  on  the  body  (for  j  |  of  orde'"  unity)  or  cm  the 

b 

surface  interactions  of  the  charged  particles  with  the  body.  The  solu¬ 
tions  in  this  Memorandum  for  jz|  ~  1  then  can  be  applied  any  body 
of  circular  cross  section. 

If  the  shape  of  the  body  in  the  ’■'lane  normal  to  V  is  other  than 

Q 

circular,  we  can  modify  our  solutions  for  Z  >  1. 

All  the  Fourier  transforms  (Eqs.  (45)  to  (58))  can  be  modified 
by  noting  that  the  effect  of  body  shape  on  the  solution  behind  the 
body  is  contained  in  the  source  term.  If  is  the  Fourier  trans¬ 
form  of  S/ (V  -  V  }f  ,  we  have 


S 


K 


-»  — »  , 
-iK-R' 
e 


-4  . 

dR' 


(65) 


where  the  double  integral  is  over  the  body  surface  in  the  piane  normal 

to  \L.  For  out  circular  cross  section, 

B 


S 


V 


(66) 


\ 


For  other  shapes,  we  can  modify  all  the  Fourier  transforms  (Eqs.  (45)  r  / 

to  (58))  by  replacing  2tt(J  (h)/h)  by  Eq.  (65). 

i 

For  the  space  dependence  of  n^.  at  large  distances  for  arbitrary 

body  shapes,  we  modify  Eqs.  (60),  (61),  (62),  and  (64)  by  multiplying 

2 

the  expression  for  (n_  -  n  )/n  by  the  factor  A/rrRD,  where  A  is  the 

i.  on  tJ 

area  of  the  body  surface  normal  to  V  .  Thus,  for  example,  Eq.  (60)  is 

D 

given  by 
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"l 
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OC 
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OD 


_ _ _ -A _ _ _ _____ _ 

r(l  +  Ta/TJ.)(2RL  V^/Vg)|sin  (z /2Ps,  )  \Jz~  sin^  ^  +  4R^  cos^  f  sin“  (z/2R^) 

(6,) 

The  solution  within  the  maximum-rarefaction  zone  will  depend  upon 
the  shape  of  the  body.  For  any  given  geometry,  the  order  of  magnitude 

$wr7 

of  e  can  be  determined  from  Poisson’s  equation* 
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